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DEFINITION OF THE SUBJECT AND ITS 
IMPORTANCE 



Cosmic strings are linear concentrations of energy that 
form whenever phase transitions in the early universe 
break axial symmetries [l| , as originally shown by Kibble 
0, Q ■ They are the result of frustrated order in the quan- 
tum fields responsible for elementary particles and their 
interactions. For about two decades, motivation for their 
study was provided by the possibility that they could be 
behind the density inhomogeneities that led to the ob- 
served large-scale structures in the universe. Precision 
observations, particularly of the cosmic microwave back- 
ground radiation, have limited strings to a sub-dominant 
role in structure formation. Instead, the inhomogeneities 
appear to be consistent with a period of cosmological in- 
flation, but it turns out that particle-physics models of 
the early universe that predict a period of inflation very 
often also predict the generation of cosmic strings at the 
end of it [1,[5|. 

More recently, interest has been revived with the re- 
alization that there may be strong links between field 
theory cosmic strings and fundamental strings. The lat- 
ter are the supposed ultimate building blocks of mat- 
ter, and in their original context of superstring theory 
were thought to be microscopic. However, in its modern 
version — sometimes referred to as M-theory — it is pos- 
sible and perhaps even mandatory to have macroscopic 
(cosmological-sized) fundamental strings [1, 0, [1| ■ Their 
behavior is expected to be quite similar to that of field 
theory cosmic strings, although there are some impor- 
tant differences so they may in principle be observation- 
ally distinguishable. Being relics of the phase transitions 
that produced them, cosmic strings provide us with a 
unique window into the early universe. If they are sta- 
ble and survive for a significant amount of time (possibly 
even up to the present day), they may leave an imprint 
in many astrophysical and cosmological observables, and 
provide us with information on fundamental physics and 
the very early universe that would otherwise be inacces- 
sible to us. On the other hand, gaining a quantitative 
understanding of their properties, interactions, evolution 
and consequences represents a significant challenge be- 
cause of their intrinsic complexity. Their non-linearity is 
particularly noteworthy, with highly non-trivial feedback 
mechanisms between large (cosmological) and small (mi- 



croscopic) scales affecting the network dynamics. Con- 
siderable reliance, therefore, must be placed on numerical 
simulations, which are technically difficult and computa- 
tionally costly. A complementary approach is the use of 
analytic or semi-analytic models, usually to describe the 
large-scale features of the networks. 

The basic picture of the cosmological evolution of 
string networks that has emerged for the simplest (Goto- 
Nambu) networks is of a scaling solution with about 40 
long strings always stretching across each horizon vol- 
ume plus a population of loops (other string types can 
lead to a different behaviour) . It is then possible to esti- 
mate their cosmological implications quantitatively. For 
example, these strings continuously source gravitational 
perturbations on sub-horizon scales. The one parameter 
in these models is the energy scale of the phase transi- 
tion at which the strings are created. The astrophysical 
consequences of strings stem from the non-trivial gravi- 
tational field around a string Q . Particles in the vicinity 
of a static straight string feel no gravitational accelera- 
tion, because in general relativity tension is a negative 
source of gravity and, since tension equals energy per 
unit length, their effects cancel. The space-time around 
the string is locally, but not globally, flat. In fact the 
space is conical, with a deficit angle 



(1) 



where /i is the energy per unit length; the simple way 
to picture this is to imagine a plane in which an angu- 
lar wedge a has been removed and the edges glued to- 
gether. For cosmologically interesting strings the deficit 
angle ranges from a few seconds of arc to a few millionths 
of a second of arc. 



II. INTRODUCTION 

To understand the cosmological evolution and effects of 
cosmic strings we start in this section with a quick sum- 
mary of basic cosmology concepts that will be needed 
later and a description of the simplest type of cosmic 
strings in which the main features are already appar- 
ent. But, first, a warning about units: from now on we 
will set the speed of light to unity c = 1, so we mea- 
sure distances in light-travel time, and masses in units 
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of Energy (and viceversa, using E = m(?). Boltzmann's 
constant is set to unity, so temperature is measured in 
units of mass/energy (using E — KbT). Finally, we set 
Planck's constant to unity, fi — \, and measure all lengths 
and masses in units of Planck's length and mass/energy: 
Ip = 1.62 X 10^35 ^ ^ ^ 2.18 X 10-^kg = 1.22 x lO^^ 
GeV/c^. In these units, Newton's constant is given by 

The early universe is very smooth. To a very good 
approximation it is a homogeneous and isotropic space- 
time described by a single variable: the rate of expansion 
of its three-dimensional spatial sections. In Einstein's 
general relativity this spacetime is described by the flat 
Friedmann-Robertson- Walker (FRW) metric 

ds"^ = dt'^ - a^{t)[dx ■ dx] , (2) 

where x are fixed (comoving) spatial coordinates and a{t) 
is the scale factor that determines the fractional or Hubble 
expansion rate 



The time coordinate t is known as cosmological time; 
to analyse cosmic string evolution we will also need a 
different time parametrization known as conformal time, 
T. They are related by dr = dt/a{t), leading to the metric 

ds^ =a'^{T)[dT^ ^dx-dx] . (4) 

The age of the universe is currently estimated to be 
about 13.7 billion years [l3|- The universe starts very 
hot and dense and is cooled by the expansion, with the 
temperature decreasing as T{t) ~ a{t)~^. The Hubble 
expansion rate is determined by the energy contents of 
the universe. In a universe dominated by radiation or 
very relativistic matter (the hottest, earliest stages), the 
scale factor evolves as a{t) ^ t^^^ and the energy den- 
sity in radiation as Eradiation ^ ~ The en- 
ergy density of non-relativistic matter is inversely pro- 
portional to volume Pmattor ~ ^-nd eventually takes 
over (after about 4000 years), leading to a period of mat- 
ter domination, during which a{t) — t^^^ and therefore 
Pmattcr ~ ■ Morc recently -about five billion years 
ago- we have entered an epoch of accelerated expansion 
due possibly to a cosmological constant or some unknown 
form of dark energy whose energy density is constant in 
time pdark energy ~ const. Dark energy should not be 
confused with dark matter, an unknown form of matter 
whose presence we can detect through its gravitational 
effects but that does not interact with electromagnetic 
fields and so in particular does not emit light -hence the 
adjective "dark"-. In the currently accepted cosmologi- 
cal model the energy density in the universe today would 
be dominated by dark energy (about 74%), followed by 
about 22% dark matter and only about 4% of regular 
(baryonic) matter [To| . Dark matter is widely believed 
to be a particle still to be discovered. 



The universe today is far from smooth, but the struc- 
ture we observe on the scale of clusters of galaxies is con- 
sistent with the gravitational collapse of tiny primordial 
density inhomogeneities 5p/ p ^ 10~^ at the time the cos- 
mic microwave background (CMB) radiation was emit- 
ted. The CMB is the oldest radiation we observe, dating 
back to the time when the universe was only 380000 years 
old. At this epoch the primordial plasma cooled enough 
to allow the formation of the first atoms (a process known 
as recombination), and it became transparent to photons 
(which is referred to as decoupling). Before that moment, 
the photons behave like a fluid that is strongly coupled 
to the protons and electrons. An overdense region in 
the baryon fluid would like to contract but the photon 
pressure pushes it back, causing both fluids to oscillate. 
These oscillations are imprinted in the cosmic microwave 
background and can be detected today in the form of 
Doppler peaks in its power spectrum. 

The spectrum of density inhomogeneities has been ac- 
curately measured in the CMB and found to be near 
scale-invariant and of the right magnitude to produce 
the structure we observe. The perturbations to the FRW 
metric can be classified as scalar (overall changes to 
the Newtonian gravitational potential), vector (associ- 
ated with velocity and/or rotational effects) and tensor 
(transverse traceless perturbations to the spatial metric, 
such as gravitational waves). Each of these affects the 
CMB in different ways, so their relative contributions 
can in principle be observationally distinguished. Fi- 
nally, Thomson scattering of the anisotropic distribution 
of the CMB photons is particularly important during de- 
coupling and recombination, and induces a partial linear 
polarization of the scattered radiation, at a level that is 
around ten percent of the anisotropy. Detection of this 
polarization signal is at the borderline of the sensitiv- 
ity of ongoing experiments at the time of writing, but is 
expected to become standard with forthcoming experi- 
ments. 

The energy density of a network of cosmic strings in the 
linear scaling regime is ^strings ~ and therefore it re- 
mains a constant fraction of the dominant form of energy 
during matter or radiation domination. Numerical esti- 
mates for the simplest, Goto-Nambu, networks suggest 
the fraction is around IQOGp. Provided the string mass 
is not close to the Planck scale, this is small enough not 
to disturb the cosmological evolution; at the same time, 
for a broad range of values of Gp this is large enough 
to be detectable in precision experiments today. Other 
string types (see section VI) may have larger or smaller 
fractions or qualitatively different signatures. In partic- 
ular, networks that do not reach linear scaling may come 
to dominate the energy density (which rules them out) 
or to disappear completely. 

The simplest field theory model that produces cosmic 
strings has a single complex scalar field <& (this is short- 
hand for a function (^{t,x) with complex values that do 
not change under coordinate transformations). Let us 
assume that the Hamiltonian determining the field dy- 
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FIG. 1: The effective potential energy V for a simple string- 
forming field theory model. The (a) and (b) plots correspond 
to the high and low temperature configurations, respectively. 
For simplicity the complex field $ has been split into two real 
scalar fields, "l>i and $2. 



namics is invariant under an axial symmetry such as a 
phase rotation, $ $e*^. For example, take the poten- 
tial energy 
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where A is a dimensionless coupling constant and rj is an 
energy scale related to the temperature of the symmetry 
breaking transition. This has a set of degenerate ground 
states: the minimum of the potential in field space is 
the circle |$| = rj, known as the vacuum manifold. Any 
configuration $(i,a;) = const. = rje''^ with x real and 
constant is a possible ground state or vacuum, irrespec- 
tive of the value of the phase x- 

Figure [1] illustrates what happens. At high tempera- 
ture the field fluctuations are large enough to make the 
central peak around |$| = irrelevant, and the effective 
potential is symmetric and has a minimum there. As the 
temperature falls the energy will eventually be too low 
to permit fluctuations over the peak, at which point the 
field will tend to settle towards one of the ground states. 
The random choice of minimum in this condensation pro- 
cess then breaks the original axial symmetry. This is the 
case, for instance, in superfluid '*He. 

When a large system goes through a phase transition 
like this, each part of it has to make this random choice, 
which need not be the same everywhere. The minimiza- 
tion of gradient terms in the energy of the system tends to 
make it evolve towards increasingly more uniform config- 
urations, but causality (the principle that no information 
can travel faster than light) imposes that this evolution 
can only happen at a limited rate. As a result one ex- 
pects many domains, each with an uncorrelated choice of 
ground state. Where these domains meet there is some 
probability of forming linear defects — cosmic strings — 
around which the phase angle varies by 27r (or possibly 
multiples thereof). This is the Kibble mechanism. No- 



tice that the field vanishes at the string's core, so there 
is trapped potential energy (as well as gradient energy). 
These strings are known as global strings because the ax- 
ial symmetry that is broken below the phase transition 
is "global" , that is, the transformation $ $e*^ is in- 
dependent of position. 

The next step is to consider charged scalar fields in- 
teracting with an electromagnetic field. The best known 
example of a symmetry-breaking transition of this kind 
is the condensation of Cooper pairs in a superconduc- 
tor, that has the effect of making photons massive below 
the critical temperature (in this case the axial symmetry 
is of the "local" of "gauge" type). The cosmic strings 
that result are magnetic flux tubes that do not dissipate 
because the magnetic field is massive outside the string 
core. 

This type of vortex was first discussed by Abrikosov 
[m in the context of type II superconductors. Nielsen 
and Olesen JJ,] generalized these ideas to the relativistic 
quantum field theory models used in particle physics, in 
particular the Abelian Higgs model which is a relativistic 
version of the Landau-Ginzburg model of superconduc- 
tivity, governed by the action. 
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is the gauge field and $ is a complex scalar of charge q 
{q = 2e in superconductors, where $ is the Cooper pair 
wavefunction) . The second term is the usual Maxwell 
action for the electromagnetic field, = dfj,Ai, — d^A^. 
The energy per unit length of a straight, static string 
lying on the z-axis is 
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where B = d^Ay — dyA^ is the z-component of the 
magnetic field. Finite energy configurations must have 
|$| = rj (the vacuum manifold is still a circle) but the 
phase of $ is undetermined provided the gradient terms 
and the magnetic field go to zero fast enough. This condi- 
tion allows for finite energy solutions At ~ Ar = A^ = 0, 
^{r,9) ~ rye'"^, Ag{r,9) ~ n/{qr), as r — > oo, in which 
the total magnetic flux in the plane perpendicular to the 
string is quantized. 
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n is the winding number of the string. If the constants A 
and q are such that fluctuations in the scalar fleld $ and 
the gauge field have equal masses, it is possible to 
show that the string saturates an inequality of the form 

Energy per unit length > constant x jmagnetic flux| 

known as the Bogomolnyi bound jisj . In this case, paral- 
lel strings at close range exert no force on each other and 
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there are static multivortex solutions [14| . If the mass of 
the scalar excitations is lower (higher) than that of the 
gauge excitations, parallel strings will attract (repel). 

More complicated particle physics models -in partic- 
ular those describing the early universe- involve gauge 
symmetries that generalize the electromagnetic interac- 
tion, mediated by photons, to more complicated inter- 
actions such as the electroweak or Grand Unified inter- 
actions. The messenger fields that play the role of the 
photons may be massless in the early universe and be- 
come massive following a symmetry-breaking transition, 
and cosmic strings carry the magnetic flux of these other 
massive gauge fields (not the electromagnetic field). 

From a cosmological point of view, the gauge field has 
the important effect of making the gradient terms decay 
exponentially fast away from the string so the energy per 
unit length of these strings is finite. Abrikosov-Nielsen- 
Olesen strings have no long-range interactions, so their 
evolution is dominated by their tension and is well de- 
scribed in the thin string or Goto-Nambu approximation. 

Field continuity implies that a string of this kind can- 
not simply come to an end: it must form a closed loop or 
extend to infinity, and it cannot break into segments. For 
this reason, strings, once formed, are hard to eliminate. 
In the absence of energy loss mechanisms, the strings 
would eventually dominate the energy density of the uni- 
verse. On the other hand, the strings can decay into 
radiation, they may cross and exchange partners, and 
they may also cross themselves, forming a closed loop 
which may shrink and eventually disappear. The out- 
come of these competing mechanisms is that the network 
is expected to reach a scale-invariant (or scaling) regime, 
where the network's characteristic length scale is propor- 
tional to the size of the horizon. We will discuss string 
evolution in more detail in Section IV. If a random tangle 
of strings was formed in the early universe, there would 
always be some strings longer than the horizon, so a few 
would remain even today. Because cosmological phase 
transitions typically happen in the very early universe, 
cosmic strings contain a lot of trapped energy, and can 
therefore significantly perturb the matter distribution. 
To first order there is a single parameter quantifying the 
effects of strings, its energy per unit length. In the sim- 
ple relativistic strings, the mass per unit length and the 
string tension are equal, because of Lorentz invariance 
under boosts along the direction of the string (but this 
need not be true for more elaborate models, see section 
VI. Cosmic strings are exceedingly thin, but very mas- 
sive. Typically, for strings produced around the epoch of 
grand unification, the mass per unit length would be of 
order fi ~ lO^^kg m^^ and their thickness 10~^^ m. The 
gravitational effects of strings are effectively governed by 
the dimensionless parameter Gfj,, where G is Newton's 
constant. For GUT-scale strings, this is 10~®, while for 
electroweak-scale strings it is 10"^"*. 



III. STRING FORMATION 

Spontaneous symmetry breaking is a ubiquitous fea- 
ture of our theories of fundamental particle interactions. 
Cosmic strings are formed in many symmetry-breaking 
phase transitions. If the symmetry is broken from a 
group G down to a subgroup H, the set of degenerate 
vacuum or ground states is the manifold M — G/H, and 
the topology of this manifold determines the types of de- 
fect that can form. In our previous examples M was a 
circle; in general, strings can form if M contains closed 
curves that cannot be contracted within M (the techni- 
cal term is that M is not simply connected, or that its 
first homotopy group is non-trivial) . 

The Kibble mechanism described in section II relates 
the initial density of strings to the size ^ of the domains 
over which the field is correlated, 

rl 

Pstring ^ ^ ^2 

with C a constant of order one reflecting the probability 
that a strings forms when three or more domains meet. 
The correlation length cannot grow faster than the speed 
of light so in the early universe an obvious upper bound 
on ^ is the size of the horizon at the time of the phase 
transition. If the dynamics of the phase transition is 
known, ^ can be estimated more accurately. In a first 
order phase transition ^ is given by the typical distance 
between bubble nucleation sites, which depends on the 
nucleation rate. In second order phase transitions ^ de- 
pends on the critical exponents and the rate of cooling 
through the critical temperature, Tc, as shown by Zurek 

MM- 

Vortex lines or topological strings can therefore ap- 
pear in a wide range of physical contexts, from cos- 
mic strings in the early universe through disclinations in 
room-temperature nematic liquid crystals, to magnetic 
flux tubes in some superconductors and vortex lines in 
low-temperature superfluid helium. These systems pro- 
vide us with a range of opportunities to test aspects of 
the cosmic string formation and evolution scenario ex- 
perimentally. 

The Kibble mechanism in first order transitions was 
confirmed in experiments on nematic liquid crystals 
[13, Ell- The Kibble- Zurek scenario for second order 
transitions has been experimentally verified in Superfluid 
^He [l^, [l^l and in Josephson Tunneling Junction arrays 

m- 

The He-3 experiments in a rotating cryostat in 
Helsinki also confirmed the scale invariance of the ini- 
tial distribution of loops, n{R) ~ where n{R)dn is 
the number density of loops with radii between R and 
R + dR, as predicted by Vachaspati and Vilenkin [2^ . 

More recently, the formation of a defect network fol- 
lowing the annihilation of ■^He-A / '^He-B boundary lay- 
ers has been observed [23]. The precise type of defects 
is still under investigation but this system constitutes an 
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interesting analogue to the formation of strings from the 
annihilation of branes in brane inflation scenarios. 

There are also a few systems where the string dei rsity 
disagrees with the Kibble-Zurek predictions. In ^He |24| . 
the reasons are miderstood: the strings are fuzzy and 
the network does not survive long enough to be detected 
[2^ . In the case of superconducting films the results are 
somewhat inconclusive [26] and also it is not completely 
clear what the expected density of flux quanta should 
be after a temperature quench; an alternative formation 
mechanism with different vortex clustering properties has 
been proposed in [23|. In fact the formation of defects 
in systems with gauge fields is clearly very relevant to 
cosmology but is still not completely understood (see [l!] 
for a recent discussion). 



IV. STRING EVOLUTION 

The motion of a cosmic string with worldsheet coor- 
dinates cr'^ and background space-time coordinates x'^ in 
a metric g^^i/ is obtainable from a variational principle 
applied to the Goto-Nambu action [2^, [s^l 



S = /^x Area = A* / dTda\detga 



drda 



for heavy strings is only relevant during a transient pe- 
riod very early in the network's evolution). The first one 
is just a wave equation with a particular damping term 
(provided by the expansion of the universe). The damp- 
ing also has the effect of reducing the coordinate energy 
per unit a. 

As was mentioned earlier the expansion of the universe 
stretches the strings, so in the absence of energy loss 
mechanisms their energy would grow with the scale fac- 
tor and the string network would eventually become the 
dominant component of the universe's energy density — 
which would be in conflict with observational results. 
Such decay mechanisms do exist (at least for the sim- 
plest models), being ultimately due to radiation losses 
and to the fact that whenever strings interact they will 
reconnect [HI, HI]- In particular closed loops may be 
formed, and these subsequently oscillate and eventually 
decay. This decay is thought to be mainly into gravi- 
tational radiation, but other forms of radiation are also 
produced very efficiently. 

Provided the decay rate is high enough, the network 
will not have pathological consequences, but will instead 
reach a linear scaling solution, where the string density 
is a constant fraction of the background density and on 
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x^gpiBem' 
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(8) 

where /i is again the string mass per unit length, and with 
dots and primes respectively denoting derivatives with 
respect to the time-like (r) and space-like (a) coordinates 
on the world-sheet, gab is called the induced metric. We 
are interested in strings in a FRW background space-time 
(see equation (|3])) and can choose worldsheet coordinates 
that make the induced metric diagonal 
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The choice of conformal time coordinate simplifies the 
microscopic evolution equations, although as we shall see 
later on physical time is a more natural choice for the 
macroscopic evolution (see Section II for definitions of 
the two time choices). It is also useful to define the co- 
ordinate energy per unit cr, 
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Then the usual variational techniques can be used to 
show that the microscopic string equations of motion are 
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and 



2e-i^ = 0. 
a 



(11) 



(12) 



For simplicity we are neglecting effects such as cusps 
and a frictional force due to particle scattering (which 



network looks the same (in a statistical 
'ries. Scaling is in fact an attractor solution, 
as has been shown both using numerical simulations and 
analytic models. Physically, the reason for this is that 
if one has a high density of strings then the number of 
string interactions increases and therefore loop produc- 
tion becomes more efficient and the decay rate increases. 
Conversely if the density is too low then there are few 
interactions and the decay rate is correspondingly lower. 
Numerical simulations confirm this broad picture, but 
also reveal that string evolution is a complex non-linear 
process, involving non-trivial interactions between vari- 
ous different scales. 

There have been thus far two generations of numeri- 
cal simulations of Goto-Nambu cosmic string networks in 
expanding universes. The first (Albrecht and Turok [ssj . 
Bennett and Bouchet [sl], Allen and Shellard [s^) dates 
from around 1990, at the peak of the interest in cosmic 
strings as possible seeds for the large-scale structures we 
observe today. In the last few years, the renewed interest 
in strings in the context of models with extra dimensions 
led to a second generation of simulations (Martins and 
Shellard :36], Ringeval et al. |37j], Vanchurin et al. ^SSJ), 
which build upon previous knowledge and exploit the 
dramatic improvements in hardware and software in the 
intervening decade and a half to achieve a much higher 
resolution. 

A different approach is provided by full field theory 
simulations [39] . These are closer to the microphysics of 
the defects and provide unique information on the inter- 
actions of the defects and their energy loss mechanisms, 
but their shorter dynamic range means that they are not 
optimal for understanding the non-linear feedback mech- 
anisms between widely different scales which affect the 
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dynamics of the network. From this point of view they 
play a very important role as calibrators, both for Goto- 
Nambu simulations and for analytic models. One can also 
carry out Minkowski space simulations (either of Goto- 
Nambu or field theory type). Neglecting the expansion 
of the universe is numerically desirable, since such sim- 
ulations are much easier to implement and evolve much 
faster. However, the expansion plays a non-trivial role 
in the network dynamics, so these results should not be 
naively extrapolated to realistic cosmological scenarios. 

Initial conditions for the numerical simulations are usu- 
ally set up using the Vachaspati-Vilenkin algorithm . 
One often adds to this random initial velocities, since 
these tend to enhance the rate of relaxation. All simu- 
lations agree on the broad, large-scale features of string 
networks, and in particular on the fact that the linear 
scaling solution is an attractor for the evolution. In Goto- 
Nambu simulations, the initial fraction of the total en- 
ergy in the form of closed loops is around 20%, but in 
the linear scaling regime this fraction is around 50% or 
even slightly more. On the other hand, in field theory 
simulations this fraction tends to be somewhat smaller. 

The first-generation simulations suggested a dynami- 
cal picture where the long-string network lost energy to 
large, long-lived loops, with sizes of order the correla- 
tion length. Refinements had each loop self-intersecting 
into around 10 daughter loops, but loop production from 
the long strings was essentially monochromatic. The 
second-generation simulations, however, reveal a quite 
different picture. Large loops do self-intersect (and in- 
deed the number of daughter loops produced by each 
one seems to around 20), but there is also a direct pro- 
duction of large quantities of small loops from slow- 
moving long-string with fractal-like substructure. In 
other words, the loop production is in fact bi-modal. All 
three second-generation simulations agree on this broad 
picture, though not on which of the two loop production 
scales is dominant. 

The second-generation simulations present some ten- 
tative evidence for the scaling of small-scale features of 
the network. An open question is whether or not this is 
expected to happen, given that gravitational backreac- 
tion (which would provide a characteristic scale) is not 
included in any of the network simulations carried out to 
date. One possible explanation stems from the fact that 
large loops are not scaled up versions of small loops. In- 
deed, small loops tend to be nearly circular, whereas large 
loops are not only far from circular but even far from pla- 
nar. In other words, the self-intersection probability for a 
given loop depends on its size, and this may be sufficient 
to dynamically select a preferred scale. Incidentally, the 
loop fragmentation processes in these networks highlight 
the fact that there is a steady flow of energy from large 
to small scales which is entirely analogous to a Richard- 
son cascade in turbulence. (In this case energy enters 
via long strings falling inside the horizon, and leaves via 
radiative decays.) 

Figure [5] shows some relevant quantities characterizing 
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FIG. 2: Characteristic small-scale properties of cosmic string 
networks in the linear scaling regime, for matter (solid lines), 
radiation (dashed) and Minkowski spacetime (dotted) runs. 
In all plots the horizontal axis represents the logarithm of 
the physical lengthscale relative to the correlation length of 
the network. The simulations leading to these results are 
described in 36]- Top panels show the correlation function 
for the velocity vectors and the fractal dimension, bottom 
panels show coarse-grained mass per unit length and coherent 
velocity. The plotted quantities are described in the main 
text. 



Goto-Nambu string networks in the linear scaling regime. 
These are always plotted as a function of scale, relative 
to the network's correlation length. The top left panel 
shows the correlation function for the velocity vectors. A 
first striking feature is that in the expanding universe cos- 
mic string velocities are anti-correlated on scales around 
the correlation length (which is smaller than but com- 
parable to the causal horizon), but such a feature is not 
present in Minkowski space. This anti-correlation is the 
result of a 'memory' of the network for recent reconnec- 
tion events, and it is ultimately due to the damping ef- 
fect of expansion. The top right panel depicts the frac- 
tal dimension of the network: this interpolates between 
d = 1 (straight segments) on small scales and d = 2 
(Brownian network) on large scales, but it does so in a 
non-trivial way (which is again different depending on 
whether or not there is expansion) and over a wide range 
of scales. The fractal dimension evolves with time, decay- 
ing on any given physical scale: the strings continually 
become smoother on any scale, so as to minimize en- 
ergy. Finally, the bottom panels show the renormalized 
(or 'coarse-grained') string mass per unit length on the 
left, and the corresponding coherent velocity on the right 
panel — notice that the effect of expansion is to reduce the 
velocities on any given scale. 

The technical difficulty and computational cost of nu- 
merical simulation provide strong motivation for alter- 
native analytic approaches, which essentially abandon 
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the detailed statistical physics of the string network to 
concentrate on its thermodynamics. The best example 
is the velocity-dependent one-scale (VOS) model [40], 
which builds on previous work by Kibble and Bennett 
and has demonstrated quantitative success when com- 
pared with both field theory and Goto-Nambu numeri- 
cal simulations. The 'one-scale' assumption is that the 
network has a single characteristic lengthscale, which co- 
incides with the string correlation length and the string 
curvature radius. This is an approximation which can be 
tested numerically. 

The first assumption in this analysis is to localize the 
string so that we can treat it as a one-dimensional line- 
like object. This is clearly a good assumption for gauged 
strings, such as magnetic flux lines, but may seem more 
questionable for strings possessing long-range interac- 
tions, such as global strings or superfluid vortex lines. 
However, good agreement between the VOS model and 
simulations has been found in both cases. The second 
step is to average the microscopic string equations of 
motion to derive the key evolution equations for suit- 
able macroscopic quantities, specifically its energy E and 
RMS velocity v defined by 



E = ^a{T) j eda , 



J eda 



(13) 



Notice that the energy is an 'extensive' quantity but the 
RMS velocity is an averaged quantity (and the averaging 
is weighted by the coordinate energy e) . In keeping with 
the above coordinate choices, the microscopic quantities 
(on the right-hand side of both equations) are defined 
in terms of conformal time, but it turns out that the 
macroscopic evolution that we are now considering is best 
described in terms of physical time — please refer to the 
cosmology review for the explicit relation between the 
two. 

Any string network divides fairly neatly into two dis- 
tinct populations, viz. long (or 'infinite') strings and 
small closed loops. In the following we will focus on the 
long strings. The long string network is a Brownian ran- 
dom walk on large scales and can be characterised by a 
correlation length L, which can be used to replace the en- 
ergy E — pV in long strings in our averaged description, 
that is. 



(14) 



A phenomenological term must then be included to ac- 
count for the loss of energy from long strings by the pro- 
duction of loops, which are much smaller than L. A loop 
chopping ejficiency parameter c is introduced to charac- 
terise this loop production as 



to loops 



~ p 

CV — 

L 



(15) 



In this approximation, we would expect the loop param- 
eter c to be a constant; comparison with numerical sim- 
ulations suggests c ^ 0.23. 



From the microscopic string equations of motion, one 
can then average to derive the evolution equation for the 
correlation length L, 



2^ ^ 2HL(l + v^) + CV , 
dt 



(16) 



where H is the Hubble parameter defined in eq. The 
first term in is due to the stretching of the network 
by the Hubble expansion which is modulated by the red- 
shifting of the string velocity, while the second is the loop 
production term. One can also derive an evolution equa- 
tion for the long string velocity with only a little more 
than Newton's second law 



? = (1 - v') 
dt ^ ' 



k{v) 



2Hv 



(17) 



The first term is the acceleration due to the curvature of 
the strings and the second is the damping term from the 
Hubble expansion. Note that strictly speaking it is the 
curvature radius R which should appear in the denom- 
inator of the first term. In the present context we are 
identifying R = L. The function k{v) is the momentum 
parameter, defined by 



k{v) 



((l-x^)(x-u)) 



(18) 



with X the microscopic string velocity and u a unit vec- 
tor parallel to the curvature radius vector. For most rel- 
ativistic regimes relevant to cosmic strings it is sufficient 
to define it as follows: 



kr{v) 



2^/2 l-Si;*^ 



1 + 8t)6 ' 



(19) 



while in the opposite case (v — > 0), we have the non- 
relativistic limit fcg = 2\/2/t:. 

Scale-invariant attractor solutions of the form L oc t 
(or L oc H^^) together with v ~ const., only appear to 
exist when the scale factor is a power law of the form 



a(t) cx t^ , 



< /3 = const. < 1 



This condition implies that 



(20) 



(21) 



with the proportionality factors dependent on the expan- 
sion rate /3. It is useful to introduce the following pa- 
rameters to describe the relative correlation length and 
densities, defining them respectively as 



L^-ft. 



C = j-^ = ptyti. 



(22) 



By looking for stable fixed points in the VOS equations, 
we can express the actual scaling solutions in the follow- 
ing implicit form: 



k{k + c) 
4/3(1 - /3) 



fc(l-/3) 
I3{k + c) 



(23) 
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where k is the constant value of k{v) given by solving the 
second (implicit) equation for the velocity. It is easy to 
verify numerically that this solution is well-behaved and 
stable for all realistic parameter values. 

If the scale factor is not a power law, then simple scale- 
invariant solutions like ([23|) do not exist. Physically this 
happens because the network dynamics are unable to 
adapt rapidly enough to the changes in the background 
cosmology. An example of this is the transition between 
the radiation and matter-dominated eras. Indeed, since 
this relaxation to a changing expansion rate is rather 
slow, realistic cosmic string networks are strictly speak- 
ing never in scaling during the matter-dominated era. 
Another example is the onset of dark energy domina- 
tion around the present day. In this case, the network 
is gradually slowed down by the accelerated expansion, 
and asymptotically it becomes frozen in comoving coordi- 
nates. The corresponding scaling laws for the correlation 
length and velocity are L cx a and v cx a~^^. 

Despite its success in describing the large-scale features 
of string networks, the VOS model has the shortcoming of 
not being able to account for the small-scale features de- 
veloping on the strings as the network evolves, as clearly 
shown by numerical simulations. This small-scale struc- 
ture is in the form of wiggles and kinks, and can be phe- 
nomenologically characterized by its fractal properties, 
as we have sketched above. As a first analytic simplifi- 
cation, the string wiggles can be characterized through 
a renormalized string mass per unit length that is larger 
than the bare (Goto-Nambu) mass. This effectively cor- 
responds to considering a model with a non-trivial equa- 
tion of state (the relation between the string tension and 
the mass per unit length), which turns out to be one 
among a larger class of models known as elastic string 
models. This kind of description has interesting parallels 
with the coarse-graining approaches that are typical of 
condensed matter. 

A more radical approach is to explicitly abandon the 
one-scale assumption. This is done in the three-scale 
model 41|, which distinguishes between the character- 
istic lengthscale (which is simply a measure of the to- 
tal string energy in a given volume) and the persistence 
length (which is defined in terms of the invariant length 
along the string and corresponds to the correlation length 
or inter-string distance). Additionally there is a third 
lengthscale which approximately describes a typical scale 
of the small-scale wiggles. This kind of description is in 
principle highly flexible, though this can be considered a 
blessing and a curse. The downside is that one is forced to 
introduce a large number of (almost free) phcnomenolog- 
ical parameters over which one has limited control even 
when comparing the model with simulations. 

Having said that, the three-scale model does confirm, 
at least qualitatively, the expectations for the behavior 
of string networks. Scaling of the large scales (in this 
case the characteristic and persistence lengths) is found 
to be an attractor, just as in the VOS model. Depending 
on the behavior of small-scale structures, the two large 



length scales may reach scaling simultaneously or the for- 
mer may do so before the latter — a behavior that has 
been seen in numerical simulations. As for the behavior 
of the small-scale structures, their evolution timescale is 
typically slower, and generically they only reach scaling 
due to the effects of gravitational backreaction (not in- 
cluded in numerical simulations). In the absence of grav- 
itational backreaction, scaling of the small-scale char- 
acteristic length is contingent of the removal of a suf- 
ficiently large amount of small-scale structure from the 
long strings by radiation and loop production, which in 
the model is controlled by a parameter whose detailed 
behavior is not known. 

Finally, an interesting and rather different approach 
starts out with the assumption that there is a range of 
scales where stretching due to the expansion is the domi- 
nant dynamical effect, even on scales well below the cos- 
mological horizon. A sufficient condition for this is that 
one is assuming that the rate of string intercommuta- 
tions is fixed in horizon units. This turns out to be suf- 
ficient to allow the construction of a statistical-type de- 
scription based on two-point correlation functions [4^ . 
Their results are to a first approximation dependent on 
a critical exponent which physically is related to the co- 
herent string velocity on a given scale. Comparison with 
numerical simulations shows, as expected, that the best 
agreement is found around and just below the horizon 
scale. 

A second assumption is that loop production at those 
scales is sufficiently localized to be describable as a per- 
turbation. When loop production is thus folded into the 
analysis, the picture that ultimately emerges is of a com- 
plicated fragmentation cascade. In particular, this model 
provides supporting evidence for the two-population loop 
distribution picture outlined above and clearly seen in 
high- resolution simulations. There is a population of 
correlation-length sized loops, produced by direct long- 
string intercommutation, and a second population with 
sizes a few orders of magnitude below (quite possibly 
near the gravitational backreaction scale) and due to loop 
fragmentation. Whether or not the smoothing provided 
by gravitational radiation is necessary to yield scaling of 
the loop sizes is again not entirely clear at the moment, 
but it is in principle a question for which this formalism 
could provide an answer. 



V. ASTROPHYSICAL AND COSMOLOGICAL 
CONSEQUENCES 

As was mentioned in Section I, the spacetime around a 
straight cosmic string is flat. A string lying along the z- 
direction has an equation of state Pz — —p, Px = Py = 
and therefore there is no source term in the relativistic 
version of the Poisson equation for the Newtonian gravi- 
tational potential 



V^(l> = AnGip + Px+Py+Pz) = 0- 



(24) 
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FIG. 3: An illustration of the mechanism behind lensing by 
cosmic strings. The thick black dot represents a cosmic string 
perpendicular to the page. The spacetime metric around the 
string can be obtained by removing the angular wedge of 
width a and identifying the edges. An observer can thus 
see double images of objects located on a certain zone behind 
the string. This zone is enclosed by the dashed lines, while 
the solid lines depict light rays and the angular separation of 
the two images, 56, will depend on the distances of the source 
and the observer to the string as well as on the deficit angle. 
Reprinted, with permission, from [4^ . 

A straight string exhibits no analogue of the Newtonian 
pull of gravity on any surrounding matter. However, this 
does not mean the string has no gravitational impact 
at all. On the contrary, we will see that a moving string 
has dramatic effects on nearby matter or propagating mi- 
crowave background photons. It is not difficult to derive 
the spacetime metric about such a straight static string 
Q. It has the simple form 

ds^ = dt^ - dz^ - dr^ - r^de^ , (25) 

which looks like Minkowski space in cylindrical coordi- 
nates, except for the fact that the azimuthal coordinate 
has a restricted range < 6' < 2-k{1 — AGfi). That 
is, the spacetime is actually conical with a global deficit 
angle 

a = SttGh ; (26) 

where an angular wedge of width a is removed and the 
remaining edges identified. 

This deficit angle implies that the string acts as a 
cylindrical gravitational lens, creating double images of 
sources behind the string (such as distant galaxies), with 
a typical angular separation 60 of order a and no distor- 
tion ^] . This is illustrated in Figure [31 A long string 
would yield a distinctive lensing pattern. We should ex- 
pect to see an approximately linear array of lensed pairs, 
each separated in the transverse direction. In each lens- 
ing event the two images would be identical and have 
essentially the same magnitude. (Except if we happen 



to see only part of one of the images.) This is a very 
unusual signature, because most ordinary gravitational 
lenses produce an odd number of images of substantially 
different magnitudes. A number of string lensing event 
candidates have been discussed in the past, but no con- 
firmed one is currently known. 

However, the above simple picture is complicated in 
practice by the fact that cosmic strings are not gener- 
ally either straight or static. Whenever strings exchange 
partners kinks are created that straighten out only very 
slowly, so we expect a lot of small-scale structure on the 
strings. Viewed from a large scale, the effective tension 
and energy per unit length will no longer be equal. Since 
the total length of a wiggly string between two points 
is greater, it will have a larger effective energy per unit 
length, U, while the effective tension T, the average lon- 
gitudinal component of the tension force, is reduced, so 
T < fi < U . This means that there is a non-zero gravi- 
tational acceleration towards the string, proportional to 
U — T. Moreover, the strings acquire large velocities, gen- 
erally a significant fraction of the speed of light, which 
introduces further corrections to the deficit angle. 

Another effect is the formation of over-dense wakes be- 
hind a moving cosmic string [i^. When a string passes 
between two objects, these are accelerated towards each 
other to a velocity 

ui_ = AttGuv , (27) 

where v is the string velocity. Matter therefore collides 
in a sheet-like structure, leaving a wake behind the mov- 
ing string. This was the basic mechanism underlying the 
formation of large-scale structures in cosmic string mod- 
els. This model has significant attractions, such as the 
early formation of nonlinear structures, and one can get 
a good match to the observed galaxy power spectrum in 
models with a large cosmological constant. However, as 
we shall discuss, it fails to reproduce the power spectrum 
of CMB anisotropies observed by COBE, WMAP and 
other experiments; cosmic strings, therefore, can only 
play a subdominant role in structure formation (albeit 
still significant, at the ten to twenty percent level). Cos- 
mic strings create line-like discontinuities in the cosmic 
microwave background signal [H, \^ . For the same rea- 
son that wakes form behind a cosmic string, the CMB 
source on the surface of last scattering is boosted towards 
the observer, so there is a relative CMB temperature shift 
across a moving string (a red-shift of the radiation ahead 
of it, and a blue-shift of that behind), given by 

ST 

— - 8ttGhv± . (28) 

where v± is the component of the string velocity normal 
to the plane containing the string and the line of sight. 
This is known as the Kaiser-Stebbins effect. This simple 
picture is again complicated in an expanding universe 
with a wiggly string network and relativistic matter and 
radiation components. The energy- momentum tensor of 
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FIG. 4: The CMB temperature power spectrum contribution 
from cosmic strings, normalized to match the WMAP data at 
£ = 10, as well as the best-fit cases from inflation only (model 
PL) and inflation plus strings (PL+S). These are compared 
to the WMAP and BOOMERANG data. The lower plot is 
a repeat but with the best-fit inflation case subtracted, high- 
lighting the deviations between the predictions and the data. 
Reprinted, with permission, from [4y]. 

the string acts as a source for the metric fluctuations, 
which in turn create the temperature anisotropics. The 
problem can be recast using Green's (or transfer) func- 
tions which project forward the contributions of strings 
at early times to today [4^ . The actual quantitative so- 
lution of this problem entails a sophisticated formalism 
to solve the Boltzmann equation and then to follow pho- 
ton propagation along the observer's line of sight. At 
the time of writing, the most recent comparisons [1^ be- 
tween full-sky maps of cosmic string-induced anisotropics 
and WMAP data yield a cosmological constraint on the 
models with 

G^i < few X IQ-'^ , (29) 

with only a weak dependence on the background 
cosmology — in particular, on the magnitude of the cos- 
mological constant. 

Apart from their scale-invariance (which follows from 
the network's attractor scaling solution discussed in 
the previous section), cosmic defect- induced fluctuations 
appear to be the antithesis of the standard inflation 
paradigm, because they are causal or active (they are 
generated inside the horizon and over an extended pe- 
riod), there are also large vector and tensor contributions, 
and they are distinctly non-Gaussian. All these charac- 
teristics leave clear signatures in the cosmic microwave 
background angular power spectrum, chief of which is a 
much broader primary Doppler peak and little evidence 
of secondary oscillations. Unlike inflation, active de- 
fect sources act incoherently with extra large-scale power 
from the vectors and tensors. Moreover, their isocurva- 



ture nature provides a partial explanation for why the 
broad primary peak ends up at larger multipoles (typi- 
cally £ ^ 300 as opposed to ^ ~ 200 in a flat cosmology. 
Isocurvature perturbations correspond to equal and op- 
posite perturbations in the radiation and matter densi- 
ties -as opposed to equal fractional perturbations in the 
number densities of the two components for adiabatic 
perturbations-. On the basis of knowledge from present 
simulations, therefore, cosmic defects alone are extremely 
unlikely to have been the seeds for large-scale structure 
formation. However, they cannot be ruled out entirely. 
For example, admixtures of inflationary power spectra 
with significant cosmic defect contributions (at a level 
around 20%, see Figured]) do provide a satisfactory fit 
to present data. This is interesting among other reasons 
because it is the sort of level at which the non-Gaussian 
signatures of cosmic strings should still be discernible, 
although their distinct line-like discontinuities are only 
clearly identifiable on small angular scales around a few 
arc minutes. 

Accelerated cosmic strings are sources of gravitational 
radiation [soj . Consequently, a network of long strings 
and closed loops produces a stochastic gravitational wave 
background [5l| over a wide range of frequencies and with 
a spectrum which (at least to a first approximation) has 
equal power on all logarithmic frequency bins. Another 
distinctive signal would come from the cusps, the points 
at which the string instantaneously doubles back on it- 
self, approaching the speed of light. Such an event gen- 
erates an intense pulse of gravitational and other types 
of radiation, strongly beamed in the direction of motion 
of the cusp [13] ■ If massive cosmic strings do indeed 
exist, both these pulses and the stochastic background 
are likely to be among the most prominent signals seen 
by the gravitational-wave detectors now in operation or 
planned, in particular LIGO and LISA. 

A stringent, though indirect, limit on the string energy 
per unit length comes from observations of the timing of 
millisecond pulsars. Gravitational waves between us and 
a pulsar would distort the intervening space-time, and so 
cause random fluctuations in the pulsar timing. The fact 
that pulsar timing is extremely regular places an upper 
limit on the energy density in gravitational waves, and 
hence on the string scale. The upper limit |53j is of order 
Gfi < 10^^, though there is still considerable uncertainty 
because this depends on assumptions about the evolution 
of small-scale structure. 

Although gravitational waves are thought to be the 
main decay byproduct of the evolution of the simplest 
cosmic string networks, direct decay into particle radi- 
ation is extremely efficient and there are claims that it 
could be the dominant energy-loss mechanism responsi- 
ble for scaling [39]. In more complicated models, there 
are certainly other decay channels. If the strings are 
global (rather than local), then they will preferentially 
produce Goldstone bosons instead. In axion models, 
these Goldstone bosons acquire a small mass and become 
the axions (a prime dark matter candidate). One can es- 
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timate the number density of axions using analytic mod- 
els of cosmic string evolution such as the VOS model. In 
another class of string models known as superconducting 
(since they have additional current-like degrees of free- 
dom), then the decay products can include electromag- 
netic radiation. 

Finally, we should mention the claims that cosmic 
strings could be responsible for a number of high energy 
astrophysical and cosmological enigmas, including ultra- 
high energy cosmic rays, gamma ray bursts, and baryo- 
genesis (the creation of the matter-antimatter asymme- 
try of the universe). Since cosmic defects can produce 
high energy particles, they could contribute to the ob- 
served cosmic ray spectrum, notably at ultra-high ener- 
gies E > lO^^GeV where the usual acceleration mech- 
anisms seem inadequate [sj]. Many ideas have been 
explored, such as particle emission from cosmic string 
cusps, but most have been found to produce a particle 
flux well below current observational limits. However, 
among the interesting scenarios deserving further study 
are those with hybrid defects (such as monopoles con- 
nected by strings) or vortons. 



VI. FIELD THEORY STRINGS WITH MORE 
DEGREES OF FREEDOM 

Kibble's original idea was to consider strings in grand 
unification scenarios, in which the strong and electroweak 
forces become unified at an energy scale of around 
2Qi5-i6QgY^ More recent studies have shown that prac- 
tically any viable Supersymmetric Grand Unified Theory 
has a pattern of symmetry breaking transitions that leads 
to the possibility of cosmic string formation at some point 
in its history 5]. In particular, in these models the infia- 
tionary phase generically ends with a phase transition at 
which strings are produced. More recent studies suggest 
that this is also a feature of brane inflation models. 

So far we have discussed the fleld theory realization 
of the simplest model of cosmic string, the Abrikosov- 
Nielsen-Olesen string, in which the mass per unit length 
equals the tension and there is no internal structure apart 
from the magnetic field. In fact, the situation can be 
much more complex in the early Universe, and realistic 
particle physics models lead to networks with a much 
richer phenomenology. The added complexity makes 
these networks much harder to study, whether by ana- 
lytic or numerical methods, and consequently they are 
not as well understood as the simplest case. 

We can only give here a brief description of the pos- 
sible complications. The list below is not complete and 
furthermore there are strings that fit more than one cat- 
egory. For a more detailed discussion we refer the reader 
to the reviews by Vilenkin & Shellard [H, Hindmarsh & 
Kibble [HI], Carter (sH], and Achiicarro and Vachaspati 
[stI ]. where references to the original literature can also 
be found. 



A. Wiggly strings; varying tension string networks; 
cycloops 



The name wiggly strings is sometimes used to refer to 
any type of string whose mass per unit length is different 
from its tension. We already mentioned in Section V 
that small structure (wiggles) on the string produces a 
renormalized effective mass per unit length [/ > /i and 
an effective tension T < fi. There are other effects that 
can affect the mass and tension, for instance the presence 
of currents along the strings. 

A particular kind of small structure is found in ex- 
tra dimensional models. If spacetime has more than the 
three spatial dimensions we observe, strings may be able 
to wrap around the extra dimensions in different ways 
leading to a renormalized four-dimensional tension and 
mass per unit length. The effective tension can of course 
vary along the strings. In extreme cases, the extra di- 
mensional wrapping effects concentrate around certain 
points along the string which behave like 'beads' (see hy- 
brid networks below) and are called cycloops. 

As discussed towards the end of Section IV, the addi- 
tional degrees of freedom (which can be thought of as a 
mass current) make the evolution of the networks highly 
non-trivial. The one-scale assumption is no longer jus- 
tified: the correlation length, inter-string distance and 
string curvature radius become distinct lengthscales. De- 
pending on the exact interplay between the bare strings 
and the mass current wiggles, these lengthscales can 
evolve differently, and some of them might be scaling 
while the others are not. The presence of extra dimen- 
sions provides a further energy flux mechanism (as energy 
may be lost into or gained from the extra dimensions) 
which will affect the string dynamics, but at the time of 
writing its exact effects have not been studied in detail. 



B. Non— topological / embedded / electroweak / 
semilocal strings 

In the Abrikosov-Nielsen-Olesen case, the scalar field 
is zero at the core of the string, and the symmetry 
is unbroken there. The zero field is protected by the 
topological properties of the vacuum manifold (the non- 
contractible circle) and the string is called topological. In 
more realistic models, the criterion for topological string 
production is a non-simply connected vacuum manifold, 
however complicated. These strings are unbreakable and 
stable. 

On the other hand there are examples in which there 
is no topological protection but the strings are never- 
theless stable. The scalar field configuration at the core 
can be deformed continuously into a ground state, so 
these non-topological strings can break, their magnetic 
flux can spread out, or be converted to a different type 
of flux. But whether this happens is a dynamical ques- 
tion that depends on the detailed masses and couplings 
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of the particles present, on the temperature, etc. 

The best studied examples of non-topological strings 
look like Abrikosov-Nielsen-Olesen strings 'embedded' 
in a larger model such as the Glashow-Salam- Weinberg 
model of electroweak interactions jSSl . [59| . These elec- 
troweak strings carry magnetic flux of the Z boson, but 
the strings would only be stable for unphysical values of 
the Z boson mass. They are closely related to semilocal 
strings, another example of embedded strings where the 
symmetry breaking involves both local and global sym- 
metries intertwined in a particular way. For low scalar 
mass these can be remarkably stable. 

In general, non-topological strings are not resilient 
enough for the networks to survive cosmological evolu- 
tion. If the strings are unstable to spreading their mag- 
netic flux, the network will not form. If the strings 
are breakable the network may form initially but it will 
quickly disappear (see hybrid networks below for a con- 
crete example). A remarkable exception to this rule are 
semilocal strings with very low scalar mass: the network 
forms as a collection of segments which then grow and 
reconnect to form longer strings or loops. These evolve 
like a network of Abrikosov-Nielsen-Olesen strings plus a 
small population of segments and there is some evidence 
of scaling [gOl • 

C. Dressed / superconducting strings / vortons 

In realistic particle physics models, a stable string will 
trap in its core any particles or excitations whose mass 
is lower inside due to the interactions with the scalar 
field. These dressed strings have a more complicated core 
structure. In extreme cases, the mass of these trapped 
particles is zero in the core and they lead to persistent 
currents along the strings, which are then known as su- 
perconducting 61]. In some cases, the decay of a loop of 
superconducting string can be stopped by these currents, 
leading to long-lived remnants called vortons [6^ that 
destroy scaling; typical vortons will either dominate the 
energy density of the universe (contrary to observations) 
or contribute to the dark matter if they are sufficiently 
light. 

D. Hybrid networks 

Hybrid networks contain more than one type of defect, 
such as for instance strings of different kinds or compos- 
ite defects combining strings, monopoles and/or domain 
walls. 



1. Composite defects: 

The production of strings may be accompanied by the 
production of other defects such as monopoles or domain 
walls, before or afterwards, that change the behaviour of 



the network as a whole. These networks can have rad- 
ically different scaling properties — in particular, linear 
scaling may not exist at all. Consider for instance a se- 
quence of breakings of the form G ^ H ^ K in which 
a symmetry group G first breaks down to a subgroup H 
which subsequently breaks to an even smaller subgroup 
K at a lower temperature. Two cases are particularly 
relevant for strings: 

The first breaking produces stable magnetic 
monopoles, the second confines -totally or partly- 
the magnetic field to flux tubes (strings) leading to a 
network of monopoles connected by strings. This can 
happen either as string segments, with monopoles at 
the ends, which eventually contract and disappear or 
as a network of strings carrying heavy "beads" (the 
monopoles) which can lead to a scaling solution. 

The first breaking produces stable strings, the second 
makes domain walls attached to the strings (e.g. in axion 
models). The network is made of pancake-like structures 
that contract under the wall tension and eventually dis- 
appear, although in some cases there may be long-lived 
remnants. 



2. Non-abelian / (p,q) strings 

Another type of hybrid network contains different 
types of strings whose intercommutation leads to three- 
point junctions and bridges. These networks are also very 
different from the simplest ones but the current consensus 
is that they also seem to reach a scaling solution during 
cosmological evolution. 

In the non-abelian case, the magnetic flux carried by 
the string is not just a number but can have different 
internal "orientations" . These become relevant when the 
strings cross, limiting the ways in which they can recon- 
nect. 

Hybrid networks containing several interacting string 
types are also found in superstring models (see next sec- 
tion). The most interesting type, usually referred to as 
(p, q) strings, contains two types of string each carrying 
different type of flux that is separately conserved: funda- 
mental and solitonic or D-strings, roughly corresponding 
to electric and magnetic flux tubes. The numbers p and 
q refer to the units of each kind of flux carried by the 
strings. Since the mass per unit length depends on these 
fluxes, (p, q) networks are expected to have a hierarchy 
of different tensions, as well as junctions and bridges. In 
fact, junctions and bridges will also form in any model 
in which parallel strings have an attractive interaction, 
such as Abrikosov-Nielsen-Olesen strings with extremely 
low scalar to vector mass ratios. 

The existence of string junctions and the hierarchy of 
string tensions make the evolution of these networks con- 
siderably more complicated than that of the simple Goto- 
Nambu strings. Relatively simple analyses suggest that 
the heavier strings with gradually decay into the lighter 
ones, and scaling is eventually reached for the strings at 
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the low end of the spectrum (the heavier ones eventually 
disappear), although this is still under discussion. Naive 
expectations that the network might be slowed down to 
non-relativistic speeds and eventually freeze have so far 
not been supported by the (admittedly simplistic) numer- 
ical simulations performed so far. Further work is needed 
to understand the general conditions under which scaling 
is (or is not) an attractor. 



VII. COSMIC SUPERSTRINGS 

Superstring theory is to date the only candidate model 
for a consistent quantum theory of gravity that includes 
all other known interactions. In string theory, the funda- 
mental constituents of nature are not point-like particles 
but one-dimensional "strings" whose vibrational modes 
produce all elementary particles and their interactions. 
Two important features of the theory are supersymme- 
try (a symmetry between bosonic and fermionic excita- 
tions that keeps quantum effects under control) and the 
presence of extra dimensions above the four spacetime 
dimensions that we observe. 

It is not yet known how to formulate the theory 
in its full generality but some weak-coupling regimes 
are well understood. In these, the fundamental strings 
live in a 10-dimensional spacetime, of which 6 dimen- 
sions are "compactified" , resulting in an effective 4- 
dimensional spacetime we live in. There is another 
regime, M-theory, in which the fundamental objects 
are two-dimensional "membranes" and the background 
spacetime is 11-dimensional. These regimes are related to 
one another by duality transformations that interchange 
the role of fluctuation quanta and non-perturbative, 
soliton-like states (branes), so the expectation is that all 
regimes are different limits of a unique, underlying the- 
ory usually referred to as superstring/M-theory, or just 
M-theory for short. 

Before the discovery of D-branes, the "solitons" of su- 
perstring theory, the question of whether fundamental 
superstrings could ever reach cosmological sizes was anal- 
ysed and the possibility discarded [6^. There were two 
main problems. First, the natural mass per unit length 
of fundamental strings is close to the Planck scale and 
would correspond to deficit angles of order 27r, which 
would have been observed. Second, the strings were in- 
herently unstable to either breaking or -depending on the 
type of string- becoming the boundary of domain walls 
that would quickly contract and disappear. The discov- 
ery of branes and their role in more exotic compacti- 
fications where the six compact dimensions have strong 
gravitational potentials (and redshifts) have changed this 
picture. It is now believed that networks of cosmic su- 
perstrings could be a natural outcome of brane-antibrane 
annihilation, especially if the branes are responsible for 
a period of cosmic inflation 0, IH, [H, . 

An important difference with previous scenarios is that 
these strings are located in regions of the compactified 



dimensions with very strong gravitational redshift ef- 
fects ( "warping" ) that reduce the effective mass per unit 
length of the strings to a level with deficit angles in the 
region of 10~^^ to 10~^, compatible with current observa- 
tions. Another important difference is a much lower prob- 
ability that the strings intercommute when they cross, 
estimated to be 10"'^ to 10^^, depending on the type 
of strings. The lower intercommutation rates lead to 
much denser networks. Estimates of the corresponding 
enhancement in the emission of gravitational radiation 
by cusps puts these strings in a potentially observable 
window by future gravitational wave detectors [13, HI] . 

The networks are hybrid, consisting of fundamen- 
tal strings and D-strings, the latter being either one- 
dimensional D-branes or perhaps the result of a higher 
dimensional D-brane where all but one dimension are 
wrapped around some "holes" (cycles) in the compacti- 
fied space. There may also be cycloops. 

As in the case of hybrid field theory strings, whether 
or not superstring networks eventually reach a scaling 
regime is an open question. Analytic studies and numer- 
ical simulations of simplified cases suggest that scaling is 
certainly possible, though contingent on model parame- 
ters that at the time of writing are not well understood. 
In this case, in addition to the presence of junctions and a 
non-trivial spectrum of string tensions, a third factor can 
affect to the evolution of these networks. If the strings 
are actually higher-dimensional branes partially wrapped 
around some extra dimensions, then energy and momen- 
tum can i n p rinciple leak into or out of these extra di- 
mensions [69| . Since the effective damping force affecting 
the ordinary and extra dimensions is different, one might 
generically expect that this will be the case. Depending 
on its sign and magnitude, such an energy fiow can in 
principle prevent scaling, either by freezing the network 
(if too much energy leaks out) or by making the strings 
dominate the universe's energy density (if too much en- 
ergy leaks in, though this is less likely than the opposite 
case) . In this sense, a somewhat delicate balance may be 
needed to ensure scaling. At a phenomenological level, 
further work will be required in order to understand the 
precise conditions under which each of these scenarios oc- 
curs. At a more fundamental level, it is quite likely that 
which of the scenarios is realized will depend on the un- 
derlying compactifications and/or brane inflation models, 
and that may eventually be used as a discriminating test 
between string theory realizations. 



VIII. FUTURE DIRECTIONS 

One of the most exciting prospects is the discovery 
of magnetic-type CMB polarization (usually referred to 
as B-modes) as this would reveal the presence of vector 
and/or tensor modes. Cosmic string models may be fur- 
ther constrained in the near future because B-modes are 
predicted to have amplitudes comparable to the electric- 
type E-modes (at large angular scales). At high reso- 
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lution, one could also hope to observe defects directly 
through the B-mode signal, against a relatively unper- 
turbed background. Conversely, the detection of vec- 
tor modes would provide strong evidence against infla- 
tion without cosmic defects. Polarization data will also 
strongly constrain a significant isocurvature contribution 
to the mainly adiabatic density fluctuations. Isocurva- 
ture perturbations can be a signature of more compli- 
cated physics during inflation, such as the effects of two 
or more scalar fields, or the formation of defects at the 
end of inflation. 

Ongoing and future CMB experiments, especially at 
high resolution, will be probing the degree of Gaussian- 
ity of the primordial fluctuations. The detection of sig- 
nificant and unambiguous non-Gaussianity in the pri- 
mary CMB signal would be inconsistent with simple (so 
called single field slow-roll) inflation. More general infla- 
tionary models can accommodate certain types of non- 
Gaussianity, and one can also envisage non-Gaussianity 
from excited initial states for inflation. It is interesting 
to note that given the existing bounds on Gfi, current 
CMB experiments do not have the sensitivity or resolu- 
tion to detect cosmic string signatures directly, in partic- 
ular the Kaiser-Stebbins effect in CMB maps. However, 
with high-resolution sensitivities becoming available in 
the near future, direct constraints (or detections) will 
be possible. This is just one example of the interest- 
ing new science that future high-resolution CMB exper- 
iments might uncover in the years ahead. In particular, 
ESA'a Planck Surveyor [See www.rssd.esa.int/Planck/], 
scheduled for launch in late 2008, may be able to provide 
significant breakthroughs. 

A deeper understanding of the evolution and conse- 
quences of string networks will require progress on both 
numerical simulations and analytic modellings. At the 
time of writing there is still no numerical code that 
includes all the relevant physics, even for the simplest 
(Goto-Nambu) strings. Inclusion of gravitational back- 
reaction is particularly subtle, and may require com- 
pletely new approaches. The expected improvements in 
the available hardware and software will allow for simu- 
lations with much longer evolution timespan and spatial 
resolution, which are needed in order to understand the 
non-linear interactions between large and small scales all 
the way down to the level of the constituent quantum 
fields. This in turn will be a valuable input for more de- 
tailed analytic modelling, that must accurately describe 



the non-trivial small-scale properties of the string net- 
works as well as the detailed features of the loop popula- 
tions. Better modelling is also needed to describe more 
general networks — three crucial mechanisms for which at 
present there is only a fairly simplistic description are the 
presence of junctions, a non-trivial spectrum of string 
tensions, and the flow of energy-momentum into extra 
dimensions. 

At a more fundamental level, a better understanding 
of the energy loss mechanisms and their roles in the evo- 
lution of the networks is still missing JXi] and it will re- 
quire new developments in the theory of quantum fields 
out of equilibrium. Such theoretical developments are 
also needed to understand defect formation in systems 
with gauge fields, and could be tested experimentally in 
superconductors. 

The early universe is a unique laboratory, where the 
fundamental building blocks of nature can be probed un- 
der the most extreme conditions, that would otherwise be 
beyond the reach of any human-made laboratory. Cos- 
mic strings are particularly interesting for this endeavour: 
they are effectively living fossils of earlier cosmological 
phases, where physical conditions may have been com- 
pletely different. The serendipitous discovery of cosmic 
defects or other exotic phenomena in forthcoming cosmo- 
logical surveys would have profound implications for our 
understanding of cosmological evolution and of the phys- 
ical processes that drove it. The search continues while, 
in the meantime, the absence of cosmic string signatures 
will remain a powerful theoretical tool to discriminate be- 
tween fundamental theories. The possibility that some- 
thing as fundamental as superstring theory may one day 
be validated in the sky, using tools as mundane as spec- 
troscopy or photometry, is an opportunity than neither 
astrophysicists nor particle physicists can afford to miss. 
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